Dark energy and dark matter as curvature effects 
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Astrophysical observations are pointing out huge amounts of "dark matter" and "dark energy" 
needed to explain the observed large scale structures and cosmic accelerating expansion. Up to 
now, no experimental evidence has been found, at fundamental level, to explain such mysterious 
components. The problem could be completely reversed considering dark matter and dark energy 
as "shortcomings" of General Relativity and claiming for the "correct" theory of gravity as that 
derived by matching the largest number of observational data. As a result, accelerating behavior of 
cosmic fluid and rotation curves of spiral galaxies are reproduced by means of "curvature effects". 
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The impressive amount of unprecedented quality data 
of last decade has made it possible to shed new light 
on the knowledge of the Universe. Type la Supernovae 
(SNela), anisotropics in the cosmic microwave back- 
ground radiation (CMBR), and matter power spectrum 
inferred from large galaxy surveys represent the strongest 
evidences for a radical revision of the Cosmological Stan- 
dard Model. In particular, the concordance ACDM model 
predicts that baryons contribute only for ~ 4% of the 
total matter - energy budget, while the exotic cold dark 
matter (CDM) represents the bulk of the matter con- 
tent 25%) and the cosmological constant A plays the 
role of the so called dark energy (~ 70%). Although 
being the best fit to a wide range of data, the ACDM 
model is severely affected by strong theoretical short- 
comings that have motivated the search for alternative 
approaches. Dark energy models mainly rely on the im- 
plicit assumption that Einstein's General Relativity is in- 
deed the correct theory of gravity. Nevertheless, it is con- 
ceivable that both cosmic speed up and dark matter rep- 
resent signals of a breakdown in our understanding of the 
laws of gravitation so that one should consider the pos- 
sibility that the Hilbert - Einstein Lagrangian, linear in 
the Ricci scalar R, should be generalized. Following this 
approach, the choice of the effective gravity Lagrangian 
can be derived by means of the data and the " economic" 
requirement that no exotic ingredients have to be added. 
This is the underivingphilosophy of what is referred to 
as f{R) gravity [IlSlESPl- From a theoretical stand- 
point, several issues from fundamental physics (quantum 
field theory on curved spacetimes, M-theory etc.) sug- 
gest that higher order terms must necessarily enter the 
gravity Lagrangian. On the other side, Solar System ex- 
periments show the validity of Einstein's theory at these 
scales. In other words, the PPN limit of /(i?)-models 
must not violate the experimental constraints on Edding- 
ton parameters. A positive answer to this request has 
been recently achieved for several f{R) theories 0, nev- 
ertheless it has to be remarked that this debate is far to 
be definitively concluded. Although higher order grav- 
ity theories have received much attention in cosmology. 



since they are naturally able to give rise to accelerating 
expansions (both in the late and the early universe |1,0]), 
it is possible to demonstrate that f{R) theories can also 
play a major role at astrophysical scales. In fact, modi- 
fying the gravity Lagrangian can affect the gravitational 
potential in the low energy limit provided that the 
modified potential reduces to the Newtonian one on the 
Solar System scale. In fact, a corrected gravitational po- 
tential could offer the possibility to fit gala xy rotation 
curves without the need of dark matter jlOl . In ad- 
dition, one could work out a formal analogy between the 
corrections to the Newtonian potential and the usually 
adopted dark matter models. The choice of an analytic 
function in term of Ricci scalar is physically corroborated 
by the Ostrogradski theorem 13], which states that this 
kind of Lagrangian is the only viable one which can be 
considered among the several that can be constructed 
by means of curvature tensor and possibly its covariant 
derivatives. The field equations of this approach, recast 
in the Einstein form, read : 

= - MR = + T'JplfiR) (1) 



where the prime denotes derivative with respect to i?, 

M 
1/3 



the standard matter stress - energy tensor and 



(2) 



defines the curvature stress - energy tensor. The presence 
of the terms f'{R)-^i, renders the equations of fourth or- 
der, while, for f{R) = R, Eqs.lQJ reduce to the stan- 
dard second -order Einstein field equations. As it is 
clear from Eq.(Q, the curvature stress - energy tensor 
formally plays the role of a further source term in the 
field equations which effect is the same as that of an 
effective fluid of purely geometrical origin. Depending 
on the scales, it is such a curvature fluid which can 
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play the role of dark matter and dark energy. From 
the cosmological viewpoint, in the standard framework 
of a spatially flat homogenous and isotropic Universe, 
the cosmological dynamics is determined by its energy 
budget through the Friedmann equations. In particu- 
lar, the cosmic acceleration is achieved when the r.h.s. 
of the acceleration equation remains positive. In phys- 
ical units, we have aja = — (l/6)(ptot + "^Vtot) ^ where 
a is the scale factor, H = a/a the Hubble parameter, 
the dot denotes derivative with respect to cosmic time, 
and the subscript tot denotes the sum of the curvature 
fluid and the matter contribution to the energy density 
and pressure. From the above relation, the accelera- 
tion condition, for a dust dominated model, leads to: 



where 
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Ptot 



so that 



a key role is played by the effective quantities : 



f'(.R) 



ifiR) - Rf'iR)] - 3HRf"{R) 



(3) 



Rf'iR) + R Rf"'{R) - Hf"{R) 
+ - L_ ., 1. (4) 



[f{R)-Rr{R)]/2-3H Rf'iR) 



As a direct simplest choice, one may assume a power - 
law form /(i?) = /o^", with n a real number, which 
represents a straightforward generalization of the Ein- 
stein General Relativity in the limit n = 1. One can find 
power -law solutions for a(i) providing a satisfactory fit 
to the SNela data and a good agreement with the esti- 
mated age of the Universe in the range 1.366 < n < 1.376 
[T^. It is worth noting, that even an inverse approach 
for the choice of /(i?) is in order. Cosmological equations 
derived from ^ can be reduced to a linear third order 
differential equation for the function /(i?(z)), where z is 
the redshift. The Hubble parameter inferred from 

the data and the relation between z and R can be used to 
finally work out /(i?) 15j. In addition, one may consider 
the expression for Hiz) in a given dark energy model as 
the input for the above reconstruction of /(i?) and thus 
work out a /(i?) theory giving rise to the same dynamics 
as the input model. This suggests the intriguing possi- 
bility of considering observationally viable dark energy 
models (such as ACDM and quintessence) only as effec- 
tive parameterizations of the curvature fluid ISj ■ 

The successful results obtained at cosmological scales 
motivates the investigation of /(i?) theories even at as- 
trophysical scales. In the low energy limit, higher order 
gravity implies a modified gravitational potential. Now, 
by considering the case of a pointlike mass m and solving 
the vacuum field equations for a Schwarzschild - like met- 
ric one gets from a theory /(i?) = foR^ the modified 
gravitational potential: 



Gm 
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which corrects the ordinary Newtonian potential by a 
power -law term. In particular, this correction sets in on 
scales larger than Tc which value depends essentially on 
the mass of the system. The corrected potential (0) re- 
duces to the standard $ (x l/r for n = 1 as it can be seen 
from the relation The generalization of Eq.® to ex- 
tended systems is straightforward. We simply divide the 
system in infinitesimal mass elements and sum up the 
potentials generated by each single element. In the con- 
tinuum limit, we replace the sum with an integral over 
the mass density of the system taking care of eventual 
symmetries of the mass distribution. Once the gravita- 
tional potential has been computed, one may evaluate 
the rotation curve v^ir) and compare it with the data. 
For the pointlike case we have : 



, Gm 

Kir) = 
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Compared with the Newtonian result — Gm/r, the 
corrected rotation curve is modified by the addition of 
the second term in the r.h.s. of Eq.(|Zj|. For < /3 < 1, 
the corrected rotation curve is higher than the Newto- 
nian one. Since measurements of spiral galaxies rotation 
curves signals a circular velocity higher than what is pre- 
dicted on the basis of the observed luminous mass and 
the Newtonian potential, the above result suggests the 
possibility that our modified gravitational potential may 
fill the gap between theory and observations without the 
need of additional dark matter. It is worth noting that 
the corrected rotation curve is asymptotically vanishing 
as in the Newtonian case, while it is usually claimed that 
observed rotation curves are fiat. Actually, observations 
do not probe Vc up to infinity, but only show that the 
rotation curve is flat within the measurement uncertain- 
ties up to the last measured point. This fact by no way 
excludes the possibility that Vc goes to zero at infinity. In 
order to observationally check the above result, we have 
considered a sample of LSB galaxies with well measured 
HI -|- Ha rotation curves extending far beyond the vis- 
ible edge of the system. LSB galaxies are known to be 
ideal candidates to test dark matter models since, be- 
cause of their high gas content, the rotation curves can 
be well measured and corrected for possible systematic 
errors by comparing 21 -cm HI line emission with op- 
tical Ha and [Nil] data. Moreover, they are supposed 
to be dark matter dominated so that fitting their rota- 
tion curves without this elusive component is a strong 
evidence in favor of any successful alternative theory 
of gravity. Our sample contains 15 LSB galaxies with 
data on both the rotation curve, the surface mass den- 
sity of the gas component and i?-band disk photometry 
extracted from a larger sample selected by de Blok & 
Bosma [la. We assume the stars are distributed in an 
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log rc 


L 


T* 


xVdof 


rms 


UGC 1230 


0.608 


-0.24 


0.26 


7.78 


3.24/8 


0.54 


UGC 1281 


0.485 


-2.46 


0.57 


0.88 


3.98/21 


0.41 


UGC 3137 


0.572 


-1.97 


0.77 


5.54 


49.4/26 


1.31 


UGC 3371 


0.588 


-1.74 


0.49 


2.44 


0.97/15 


0.23 


UGC 4173 


0.532 


-0.17 


0.49 


5.01 


0.07/10 


0.07 


UGC 4325 


0.588 


-3.04 


0.75 


0.37 


0.20/13 


0.11 


NGC 2366 


0.532 


0.99 


0.32 


6.67 


30.6/25 


1.04 


IC 2233 


0.807 


-1.68 


0.62 


1.38 


16.29/22 


0.81 


NGC 3274 


0.519 


-2.65 


0.72 


1.12 


19.62/20 


0.92 


NGC 4395 


0.578 


0.35 


0.17 


6.17 


34.81/52 


0.80 


NGC 4455 


0.775 


-2.04 


0.88 


0.29 


3.71/17 


0.43 


NGC 5023 


0.714 


-2.34 


0.61 


0.72 


13.06/30 


0.63 


DDO 185 


0.674 


-2.37 


0.90 


0.21 


6.04/5 


0.87 


DDO 189 


0.526 


-1.87 


0.69 


3.14 


0.47/8 


0.21 


UGC 10310 


0.608 


-1.61 


0.65 


1.04 


3.93/13 


0.50 



TABLE I: Best fit values of tfie model parameters from max- 
imizing the joint likelihood function £(/3, log rc, /g). We also 
report the value of T*, the /dof for the best fit parameters 
(with dof — N — 3 and A'^ the number of datapoints) and the 
root mean square arms of the fit residuals. 



infinitely thin and circularly symmetric disk with sur- 
face density = Tj,lQexp{—R/ Rd) where the central 
surface luminosity Iq and the disk scalelength Rd are ob- 
tained from fitting to the stellar photometry. The gas 
surface density has been obtained by interpolating the 
data over the range probed by HI measurements and ex- 
trapolated outside this range. 

When fitting to the theoretical rotation curve, there 
are three quantities to be determined, namely the stellar 
mass - to - light (M/L) ratio, and the theory param- 
eters (/3, Tc). It is worth stressing that, while fit results 
for different galaxies should give the same f3, Tc must 
be set on a galaxy - by - galaxy basis. However, it is ex- 
pected that galaxies having similar properties in terms 
of mass distribution have similar values of rc so that the 
scatter in Vc must reflect somewhat that on the terminal 
circular velocities. In order to match the model with the 
data, we perform a likelihood analysis determining for 
each galaxy using as fitting parameters /3, logr^ (with 
Tc in kpc) and the gas mass fraction^ fg. Considering 
the results summarized in Tabled the experimental data 
are successfully fitted by the model. In particular, for 
the best fit range of /3 {(3 — 0.58 ± 0.15), one can con- 
clude that i?" gravity with 1.34 < n < 2.41 (which well 
overlaps the above mentioned range of n interesting in 
cosmology) can be a good candidate to solve the missing 
matter problem in LSB galaxies without any dark mat- 



1 This is related to the M/L ratio as T* = [(1 - fg)Mg]/{fgLa) 
with Mg = IAMhi the gas (HI + He) mass, Aid = '^*Ld and 
Ld = 2itIqB?^ the disk total mass and luminosity. 
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FIG. 1: Best fit theoretical rotation curve superimposed to 
the data for the LSB galaxy NGC 4455 (left) and NGC 5023 
(right) . To better show the effect of the correction to the New- 
tonian gravitational potential, we report the total rotation 
curve Vc(R) (solid line), the Newtonian one (short dashed) 
and the corrected term (long dashed). 

ter. At this point, it is worth wondering whether a link 
may be found between i?" gravity and the standard ap- 
proach based on dark matter haloes since both theories 
fit equally well the same data. As a matter of fact, it 
is possible to define an effective dark matter halo by im- 
posing that its rotation curve equals the correction term 
to the Newtonian curve induced by i?" gravity. Mathe- 
matically, one can split the total rotation curve derived 
from i?" gravity as v1{r) = ^{r) + vl^^^^A''') where the 
second term is the correction one. Considering, for sim- 
plicity a spherical halo embedding an infinitely thin ex- 
ponential disk, we may also write the total rotation curve 
as vl{r) = vl^,,k{r) + vljjMir) with w^_^,,fc(r) the New- 
tonian disk rotation curve and Vc DAii''^) — GMuM(,r) /r 
the dark matter one, MoMir) being its mass distribu- 
tion. Equating the two expressions, we get : 

MdmM = 2P-Sc^<^ - moRh^Mv,f3) ■ (8) 
with 1] = r/R^i^ Sq = T^/q and: 

Jo 

with To only depending on the geometry of the system. 
Eq.© defines the mass profile of an effective spheri- 
cally symmetric dark matter halo whose ordinary rota- 
tion curve provides the part of the corrected disk rota- 
tion curve due to the addition of the curvature corrective 
term to the gravitational potential. It is evident that, 
from an observational viewpoint, there is no way to dis- 
criminate between this dark halo model and R" gravity. 
Having assumed spherical symmetry for the mass distri- 
bution, it is immediate to compute the mass density for 
the effective dark halo as p_DA/(r) = {l/ATrr^)dMDM /dr. 
The most interesting features of the density profile are 
its asymptotic behaviors that may be quantified by the 
logarithmic slope aoM — din puM /dlnr which can be 
computed only numerically as function of rj for fixed val- 
ues of P (or n). The asymptotic values at the center and 
at infinity denoted as ao and ttoo result particularly in- 
teresting. It turns out that ao almost vanishes so that in 
the innermost regions the density is approximately con- 
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stant. Indeed, ao = is the value corresponding to mod- 
els having an inner core such as the cored isothermal 
sphere and the Burkert model Moreover, it is well 
known that galactic rotation curves are typically best 
fitted by cored dark halo models therein). On the 
other hand, the outer asymptotic slope is between —3 
and —2, that are values typical of most dark halo mod- 
els in literature. In particular, for /3 = 0.58 one finds 
(ao,aoo) — (—0.002,-2.41), which are quite similar to 
the value for the Burkert model (0,-3), that has been 
empirically proposed to provide a good fit to the LSB and 
dwarf galaxies rotation curves. The values of (ao, Ofco) we 
find for our best fit effective dark halo therefore suggest a 
possible theoretical motivation for the Burkert - like mod- 
els. Now, due to the construction, the properties of the 
effective dark matter halo are closely related to the disk 
one. As such, we do expect some correlation between 
the dark halo and the disk parameters. To this aim, ex- 
ploiting the relation between the virial mass and the disk 
parameters , one can obtain a relation for the Newtonian 
virial velocity T^ir = GMyi^/ Rmr ■ 

1-13 1±£ a ^ 

We have numerically checked that Ea. (|10|l may be 
well approximated as Md oc V^^J. which has the same 
formal structure as the baryonic TuUy- Fisher (BTF) 
relation Mf, oc Vfl^^ with the total (gas -I- stars) 
baryonic mass and Vfiat the circular velocity on the 
flat part of the observed rotation curve. In order to 
test whether the BTF can be explained thanks to the 
effective dark matter halo we are proposing, we should 
look for a relation between Vmr and Vfiat- This is 
not analytically possible since the estimate of Vfiat 
depends on the peculiarities of the observed rotation 
curve such as how far it extends and the uncertainties 
on the outermost points. For given values of the disk 
parameters, we therefore simulate theoretical rotation 
curves for some values of and measure Vfiat finally 



choosing the fiducial value for Vc that gives a value 
of Vfiat as similar as possible to the measured one. 
Inserting the relation thus found between Vfiat and Vmr 
into Ea. (|10l) and averaging over different simulations, we 
finally get: logM,, = (2.88±0.04) logj^/at + (4.14±0.09) 
while observational data give [13 : log Mb = 
(2.98 ± 0.29) log Vf/at + (3.37 ± 0.13) . The slope of 
the predicted and observed BTF are in good agreement 
thus leading further support to our approach. The 
zeropoint is markedly different with the predicted one 
being significantly larger than the observed one, but 
it is worth stressing, however, that both relations fit 
the data with similar scatter. A discrepancy in the 
zeropoint may be due to our approximate treatment 
of the effective halo which does not take into account 
the gas component. Neglecting this term, we should 
increase the effective halo mass and hence Vyir which 
affects the relation with Vfiat leading to a higher than 
observed zeropoint. Indeed, the larger is Alg/M^, the 
more the point deviate from our predicted BTF thus 
confirming our hypothesis. Given this caveat, we may 
therefore conclude with confidence that i?" gravity offers 
a theoretical foundation even for the empirically found 
BTF relation. 

The results outlined along this paper are referred to a 
simple choice of f{R), while it is likely that a more com- 
plicated Lagrangian is needed to reproduce the whole 
dark sector phenomenology at all scales. Nevertheless, 
although not definitive, these achievements represent an 
intriguing matter for future more exhaustive investiga- 
tions. In particular, exploiting such models can reveal a 
useful approach to motivate a more careful search for a 
single fundamental theory of gravity able to explain the 
full cosmic dynamics with the only two ingredients we 
can directly experience, namely the background gravity 
and the baryonic matter. 

Acknowledgements We wish to thank Massimo 
Capaccioli for the useful discussions and encouragements 
which allowed to improve the paper. 



[1] S. Capozziello, Int. J. Mod. Phys. D, 11, 483 (2002) 
[2] S.M. Carroll, V. Duvvuri, M. Trodden, M. Turner, Phys. 

Rev. D, 70, 043528, 2004 
[3] E.E. Flanagan, Class. Quant. Grav., 21, 417, 2003 
[4] G. AUemandi, A. Borowiec, M. Francaviglia, Phys. Rev. 

D, 70, 103503, 2004 
[5] S. Nojiri, S.D. Odintsov, 2006, hep-th/0601213 
[6] N.D. Birrell, P.C.W. Davies, Quantum Fields in Curved 

Space, Cambridge Univ. Press, Cambridge (UK), 1982 
[7] S. Capozziello, A. Troisi, Phys. Rev. D, 72, 044022, 2005 
[8] A.A. Starobinsky, Phys. Lett. B, 91, 99, 1980 
[9] R. Kerner, Gen. Rel. Grav., 14, 453, 1982 
[10] S. Capozziello, V.F. Cardone, A. Troisi, 2006, 

astro-ph/0603522 
[11] K. Stelle, Gen. Rel. Grav. 9, 353, 1978 



[12] M. Milgrom, Astroph. Journ., 270, 365, 1983 

[13] R.P. Woodard, astro-ph/0601672 2006 

[14] S. Capozziello, V.F. Cardone, S. Carloni, A. Troisi, Int. 

J. Mod. Phys. D, 12, 1969, 2003 
[15] S. Capozziello, V.F. Cardone, A. Troisi, Phys. Rev. D, 

71, 043503, 2005; T. Multamaki and I. Vilja, Phys. Rev. 

D, 73, 024018, 2006. 
[16] W.J.G. de Blok, A. Bosma, Astron. & Astroph. 385, 816, 

2002 

[17] A. Buvkeit, Astroph. Journ., 447, L25, 1995 
[18] G. Gentile, P. Salucci, Mon. Not. Roy. Astron. Soc., 351, 
953, 2004 

[19] S.S. McCaugh, Astroph. Journ., 632, 859, 2005 



